Donnan equilibrium and the osmotic pressure of charged colloidal lattices 
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We consider a system composed of a monodisperse charge-stabilized colloidal suspension in the 
presence of monovalent salt, separated from the pure electrolyte by a semipermeable membrane, 
which allows the crossing of solvent, counterions, and salt particles, but prevents the passage of 
polyions. The colloidal suspension, that is in a crystalline phase, is considered using a spherical 
Wigner-Seitz cell. After the Donnan equilibrium is achieved, there will be a difference in pressure 
between the two sides of the membrane. Using the functional density theory, we obtained the 
expression for the osmotic pressure as a function of the concentration of added salt, the colloidal 
volume fraction, and the size and charge of the colloidal particles. The results are compared with 
the experimental measurements for ordered polystyrene lattices of two different particle sizes over 
a range of ionic strengths and colloidal volume fractions. 

PACS numbers: 82.70.Dd: 36.20.-r; 64.60.Cn 



I. INTRODUCTION 

In recent years, colloidal particles have received an 
increased attention because they constitute an interest- 
ing system from practical, experimental, and theoretical 
point of view. Large molecules immersed in a solution 
are important in various systems, from biological to in- 
dustrial. Due to the large size of the particles, a rich 
variety of experiments can be easily performed. Optical 
measurements show that some suspensions, such as the 
opals and the viruses ||], form regular lattices that 
can, in principle, exhibit melting and structural phase 
transitions The elastic rigidity of these ordered 

structures leads to unusual viscoelastic properties. The 
suspension responds to small-amplitude deformations as 
a linear viscoelastic solid, allowing for the propagation of 
the low- frequency shear waves . 

From the purely theoretical point of view, colloids are 
quite fascinating materials to study. Although thermody- 
namically identical to the usual atomic fluids and solids, 
colloids have an additional advantage in that the range 
of the interactions between the colloidal particles can 
be "manually" controlled by exploiting their interactions 
with a surrounding solvent, as well as by actually synthe- 
sizing particles which behave in a desired fashion. In gen- 
eral, the interactions between the colloids are dominated 
by a van der Waals force resulting from the quantum fluc- 
tuations of the electron charge density on the surface of a 
colloidal particle. This attractive interaction can lead to 
an aggregation of macromolecules and to their precipita- 
tion in a gravitational flcld. In many practical applica- 
tions, such as design of water-soluble paints, it is essential 
to device a mechanism which would stabilize the colloidal 
suspension against precipitation. One such mechanism is 
to synthesize colloidal particles with some acidic groups 
on their surfaces, which will be ionized upon contact with 
water. The sufficiently strong electrostatic repulsion be- 
tween the equally charged macromolecules will prevent 
the formation of clusters and stabilize the suspension 



against the precipitation. It is the aim of this paper to try 
and shed some light on the behavior of charge-stabilized 
colloids. 

When the volume fraction is not too small, the charged 
polyions form an ordered structure (bcc or fee) 0-^- 
The description of the suspension in this case becomes 
significantly more simple than that of a disordered struc- 
ture 10Hl2i, since one can take advantage of the transla- 
tional symmetry of the lattice. Thus each polyion, with 
its counterions, is inclosed in a Wigner-Seitz (WS) cell 
10-^. The thermodynamics of the system is then fully 
determined by the behavior inside one cell. In this paper 
we will present a simple theory to calculate the Donnan 
equilibrium properties of a charged colloidal lattice. The 
central problem will be to obtain the osmotic pressure 
of a colloidal crystal in equilibrium with a reservoir of 
salt. The experimental setup is described in Ref. It 
consists of an osmometer, which is made of two cells sep- 
arated by a semipermeable membrane, which allows for 
the crossing of solvent, counterions, and the microions 
of salt, but prevents the passage of polyions. This, in 
turn, leads to an establishment of a membrane potential 
and an imbalance in the number of microions on the two 
sides of the membrane, resulting in an osmotic pressure 
measured by two capillaries attached to the chambers 
containing the colloidal suspension and the pure salt so- 
lution. This osmotic pressure, measured in mm of H2O, 
will be compared to that found on the basis of our theo- 
retical calculations. 



II. OSMOTIC PRESSURE 



Since the colloidal suspension is organized in a peri- 
odic structure, it is sufflcient to consider just one iso- 
lated polyion in a salt solution inside an appropriate WS 
cell. To simplify the problem, a further approximation 
is to replace the polyhedral WS cell by a sphere of the 
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same volume V = ^R^, as represented in Fig. 1. Previ- 
ous calculations, which compare explicitly results derived 
from the spherical-cell model with those that follow from 
a periodic system of cubic symmetry, show that this ap- 
proximation provides a good description of the multibody 
system at low- volume fractions [Q. 




FIG. 1. Geometry of the spherical Wigner-Seitz cell with 
radius R. The mesoscopic polyion (with total charge —Ze 
distributed uniformly on its surface) is represented by the 
shaded sphere of radius a. The semipermeable membrane lo- 
cated at \r\ = R, represented by the dashed line, separates 
the polyelectrolyte from the pure solution of salt. The micro- 
scopic mobile ions — counterions (©), cations and anions 
(— ) — are free to move in the region |r| > a, and can cross 
the semipermeable membrane. 

Our theoretical picture, exploiting the underlying pe- 
riodicity of the colloidal suspension jl4j, consists of a 
spherical polyion of radius a with a uniform surface 
charge density. 
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residing at the center of a spherical WS cell, where Z is 
the polyion valence, e is the electronic charge and S is the 
Dirac distribution. Note that we assume that the polyion 
is negatively charged. The boundary of the WS cell cor- 
responds to the semipermeable membrane. The solvent, 
usually water, is modeled by a homogeneous continuum 
of a dielectric constant e. The mobile microions, modeled 
as point particles of charge ±e, include the counterions 
and the ions derived from the dissociation of the mono- 
valent salt, cations and anions. Inside the PFS'cell, these 
are free to move within the annulus a < \r\ < R. The 
distribution of microions is strongly determined by the 
electrostatic potential produced by the central polyion. 
In the region |r| > J? the microions are assumed to be un- 
affected by the presence of the polyion at \r\ =0 and are 
uniformly distributed . Within the functional density 
theory, the microions inside the WS cell are treated as 
an inhomogeneous fluid. At this level of approximation 



we shall not distinguish between the counterions and the 
cations of salt, so that the local density of positive par- 
ticles is p+{r). The local density of co-ions (anions) is 
p_(r). 

The effective Helmholtz free energy inside the WS cell, 
associated with the region of the colloidal suspen- 
sion, is composed of: 

(a) entropic ideal-gas contributions associated with the 
microion local densities p±{r), 
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where P = (kBT)-^ is the inverse temperature and 
A = h/{2mnkBTy^^ is the thermal de Broglie wave- 
length; 

(b) electrostatic energy terms due to the polyion- 
microions interactions. 
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(c) electrostatic energy terms due to the microions- 
microions interactions. 



F^j = e^ I (fr d^r 
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for i, j = ± ; 
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(d) and finally a (virtual) Lagrange-multiplier term, 

F'' =pD j d'r[p+{r)~p^{r)] , (5) 

which enforces the overall charge neutrality within the 
WS cell. 
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With the introduction of (^), we can treat the densities 
of microions, p±{r), as independent variables. The La- 
grange multiplier p,£) is traditionally referred to as the 
Donnan potential. After the minimization procedure 
is complete, p,D shall be eliminated using the charge- 
neutrality constraint (^. We have, then, the Helmholtz 
free energy inside the WS cell. 
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In the region outside the WS cell, the salt solution is, 
once again, charge neutral and the positive and negative 
ions are uniformly distributed. Due to this symmetry, the 
overall Coulombic contribution vanishes up to the corre- 
lation free energy, which we neglect within our treatment 
fl^ . In the absence of an external field, the densities of 



positive and negative ions within the electrolyte chamber 
are uniform and given by 



p+{\r\>R) = p^{\r\>R)=ps. 



(8) 



where ps is the (uniform) salt (cation/anion pairs) con- 
centration of the pure electrolyte in the chamber after the 
equilibrium has been established. The Helmholtz free en- 
ergy, F°^^, then reduces to that of an ideal gas, 



= Vs p+ [In (p+A3) -1]+Vs p- [In (p_ A^) - l] , 

(9) 

where Vs is the volume of the chamber containing pure 
electrolyte. 

Assuming that the region outside the WS cell corre- 
sponds to the chamber containing the pure-electrolyte 
solution, the system will achieve the Donnan equilibrium 
when the electrochemical potentials inside. 
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= ^In [p±(t')A^] ± fiD ± e^Pir) , 
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and outside the WS cell, 
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are equal. 



p^{r)=pr 



(12) 



These conditions lead to the Boltzmann distribution of 
the microion densities, 



P±{r) = PS exp [TPt^D T el3ip{r)] 



(13) 



On the other hand, the electrostatic potential ip{r) and 
the local charge densities obey the (exact) Poisson equa- 
tion, 



vV(^) = V-£;(r) 
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[qp{r) + ep+{r) -ep_(r)] , 
(14) 



where E{r) is the electric field at the point r, which is 
related to the electrostatic potential by E{r) = —Vilj{r). 
Inserting expression (^) into Eg. ([l^) we find a Poisson- 
Boltzmann-/ifce equation [|l6|j8[[7|,|5|,|l^ . The effects of the 
correlations inside the ionic atmosphere can be included 
through the use of the local or the weighted-density con- 
tribution to the free-energy Q. This will be a topic for 
the future. 

We now take advantage of the spherical symmetry of 
the system to eliminate the angular dependence of the 
equations, that is, we replace r by r |r| in Eqs. ( |l3|) 
and (|l^). The electric field also has a spherical sym- 
metry, so that E{r) = E{r)Y- Integrating the Poisson 



equation (|lj) over a sphere of radius r and using the di- 
vergence theorem, we obtain the equation for the electric 
field strength E{r), 

j (fr' V • E{r') = J dS' ■ E{r') = 4:TTr^E{r) 
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where 



a(r)= I dV [p+(r')-p-(r') 



(15) 
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Using the Boltzmann equations (|l^) for the densities of 
microions inside the WS cell, we can then write 

a{r) = Anps [exp (-Pfio) a+ir) - exp (/3/iz)) a_(r)] , 

(17) 

where a± (r) are the radial integrals associated with the 
electrostatic potential ip{r), that can be written in terms 
of the electric field strength E{r), 



<^±{r) ^ / dr'r exp[^e(3tp{r')] 
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where we have chosen the gauge in which tpir — a) — 0. 
The Lagrange multiplier pLo can now be eliminated using 
the charge-neutrality constraint over the WS cell (|), 



a{R) = Attps [exp {-(ipo) a+{R) - exp {iSpn) a-{R)] = Z , (19) 
which gives a quadratic equation for exp (— /3/i/)). Solving the associated equation and replacing in (p^), we obtain 

a+(r) 
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To perform the numerical calculations, it is convenient 
to write the equations in terms of dimensionless variables, 



f = - , Eir) = ef3\E{r) , (21) 
p± (f) = 47rA=^p± (r) , ps = AnX^ps , 
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where we have introduced the Bjerrum length. 



(22) 



charge is uniformly distributed and the pressure, up to 
the electrostatic correlations effects, is that of an ideal- 
gas, 



dF° 
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(27) 



Therefore, the osmotic pressure, defined as the difference 
between the pressures inside and outside the WS cell, is 



/3AP = /3P" 
1 



/3P° 



p+{R)+P-{R)-2ps 
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(28) 



which measures the strength of the electrostatic interac- 
tions in the colloidal suspension. 

The equation for the dimensionless electric field 
strength can be written as 
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(23) 



where a{f) is given by Eq. (20), replacing a±(r) by 
a±{f), and 87rp5 by 2ps- Together with (|20|), this is 
an integral equation which determines the reduced elec- 
tric field E{f). It is gauge invariant, since it does not 
depend on the choice for tpi^r = a), and already contains 
the boundary conditions. 



E{r = a) = 



and 



E{r = R)=^. 



(24) 



This equation can be solved numerically by replacing the 
continuous interval a < f < i? by a finite grid and it- 
erating ( p3| ) until convergence is found. In the absence 
of added salt, ps = 0, we regain the results of previous 
calculations [^]. 

The dimensionless densities of microions inside the WS 
cell are given by 



P±(^) = 



Z2-f 4p|ci+(i?)a_(i?)±Z 



X exp 



2a±(i?) 
± / df' E{f') 



(25) 



The pressure inside the cell is P™ - 

1 ■r-in j; — L 

4 ji^ dR ' after some algebra reduces to ||16[ 



/3P- = p+(i?) + p_(i?) 



(26) 



The physical interpretation of (|2^) is that at the border 
of the WS cell the electrostatic forces on the mobile ions 
vanish, since E{R) — 0, and the pressure is given by the 
ideal-gas law. On the other side of the membrane the 



III. COMPARISON WITH EXPERIMENTAL 
RESULTS 

The derivation in the previous section requires us to 
know the concentration of salt ps once equilibrium across 
the membrane has been established. Although, in princi- 
ple, this is possible to measure, it is much simpler, from 
the point of view of experiment, if the osmotic pressure 
was calculated as a function of the concentrations of salt 
present in the two chambers, and p^, before the ion 
transport is allowed to take place. Evidently it is these 
values that are directly controlled by the experimentalist. 
In this section we will present the results for the osmotic 
pressure as a function of the original concentration of salt 
added to the polyelectrolyte in the first chamber, 7i°, and 
of the initial pure-electrolyte concentration in the second 
chamber, p^, before the ion transport is allowed to take 
place. In this way our results can be directly compared 
with the measurements of Benzing and Russel . 

The initial concentration of salt in the pure-electrolyte 
chamber, p'g, will not coincide with the same quantity ps 
measured after the Donnan equilibrium is achieved. A re- 
lation between these two concentrations can, however, be 
determined through a mass-conservation condition as fol- 
lows. If and n are the (macroscopic) average salt con- 
centrations in the colloidal-suspension chamber before 
and after the Donnan equilibrium is achieved, and Vc and 
Vs are the volumes of the colloidal-suspension and the 
pure-electrolyte chambers, then the mass-conservation 
condition, 



or 



n°Vc + p%Vs ^ nVc + psVs 



n ^ n"(l + 7172) - 7ips 



(29) 



(30) 



gives a relation between the average concentrations of 
salt in the suspension after and before the equilibriimr in 
terms of the ratios 



71 = — and 
Vc 



(31) 



Since the microions are just allowed to occupy an ef- 
fective volume of the WS cell, 



(32) 



where is the volume fraction occupied by the polyion, 

(33) 



it is useful to define an effective salt concentration in the 
suspension, ng, given by 
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(34) 



This effective salt concentration is also obtained by inte- 
grating the density of negative salt particles over the WS 
cell, 



where 



1 / Z 
ns^P- = 7;[P+ + P---r7- 



p± = — \ d^rp±{r) 

a<\r\<R, 



(35) 
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Using the charge-neutrality constraint over the WS 
cell (O), we have 



P+ - P- 



Z 



(37) 



On the other hand, integrating the Boltzmann equa- 
tions (O), we obtain 



P± = ^ expiTPfiD) r dr' r'^ exp [Te/3V(r')] 
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Solving (HI) and (||), we find 
Z 
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Using now Eq. (|35|), we can relate the density of salt in 
the pure-electrolyte chamber, ps, and the effective con- 
centration of salt in the colloidal suspension, ns- After 
the equilibrium is achieved, 

(47rps)'a+(i?)«-(i?) = 4Kff + ZngV^s ■ (41) 

Substituting Eq. ( ^4|) for ns, we are left with a quadratic 
equation for ps, which can be easily solved to yield 
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(42) 



where 

1^ = nsV,B + yPsi?' = ^^R^n°{l + 7172) 



(43) 



Inserting ps into the expression for a{r), the integral 
Eq. ( psf ) can, once again, be solved iteratively to obtain 
the electric field and the density profiles for the microions 
inside the WS cell. The osmotic pressure is calculated as 
in the previous section, but now is function of the initial 
densities of salt n° and p^. 

We are now ready to compare our results for 
the osmotic pressure with the experimental measure- 
ments of Benzing and Russel for monodisperse ordered 
polystyrene lattices In their experiments Vc = Vs 
and the initial density of salt inside the colloidal suspen- 
sion was adjusted so that n'^ — p^{l — (f)) to minimize the 
ion transport. We then have 
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1, 



and 



72 



1 



(44) 



The osmotic pressure was measured in terms of the vol- 
ume fraction (j) for several initial salt concentrations n° 
for two distinct sizes of latex polyions. The experimen- 
tal results are reproduced in Figs. 2-3. From Eqs. (p8|), 
(|4^ ) and (44), we obtain the osmotic pressure for a fixed 
Z in terms of the volume fraction (j) for different initial 
concentrations of salt n° inside the colloidal suspension. 
In order to compare the results of our theoretical cal- 
culations with the measurements of Benzing and Russel 
1^, we need the valence Z of the polyions. We shall use 
Z as a fitting parameter, determined in such a way as 
to make the theoretical curve for the lowest electrolyte 
density passes through one experimental data point. For 
instance, for the latex A (a = 525 A), we chose the exper- 
imental point n'-* = 10^^ mole/dm^ and = 0.3, and we 



found for the polyion valence the value Za = 923. Using 
always the same value of Z and varying n", we obtained 
the other graphs (see Fig. 2). The same procedure was 
appHed to the latex B (a = 1830 A) — see Fig. 3 — 
where we have used the point vl^ — 10~^ mole/dm"^ and 
= 0.2 and found Zb = 2212. The general trend, ob- 
served both in the experimental measurements as well as 
in the theoretical calculations, is that the osmotic pres- 
sure was decreased by an increase in the particle size 
or the amount of added salt. Comparing our theoretical 
predictions with the experiments, one might note that al- 
though the curves for the lowest salt density provide an 
almost perfect fit to the experiment, as the electrolyte 
density is varied, such a good agreement is lost. One 
possible reason for this behavior might lie in the surface 
chemistry of the latex particles, which causes the num- 
ber of ionized sites Z to be dependent on the strength of 
electrolyte in such a way as to make Z increase with the 
added concentration of salt fl^ . 



IV. CONCLUDING REMARKS 



We have considered a system composed of a monodis- 
perse charge-stabilized colloidal suspension and monova- 
lent salt separated from the pure electrolyte by a mem- 
brane permeable to the microions and solvent, but im- 
permeable to the polyions. Introducing a spherical WS 
cell, we have studied the crystalline phase within the 
framework of the density functional theory. We have nei- 
ther considered correlations between the microscopic free 
ions, nor the polyion-polyion interaction explicitly. The 
latter was replaced by an appropriate boundary condi- 
tion, namely, the charge-neutrality constraint over the 
WS cell. 

In particular, we have determined the osmotic pres- 
sure between the two sides of the membrane as a func- 
tion of the initial concentration of added salt and the 
size, charge, and concentration of the colloidal particles. 
Using the valence of the polyion Z as a fitting param- 
eter, we compared our theoretical results with experi- 
mental measurements for ordered polystyrene lattices of 
two different sizes , over a range of ionic strengths and 
colloidal volume fractions. Our results are in agreement 
with the experimental data, which shows that the os- 
motic pressure is decreased by an increase in the parti- 
cle size or the amount of added salt, or by dilution of 
the colloidal particles in the suspension. Although the 
experimental data corresponding to the lowest salt con- 
centrations agrees reasonable well with our theoretical 
predictions, such good agreement is lost for higher elec- 
trolyte densities. However, the trend of our results agrees 
with the fact that the valence Z increases with the salt 
concentration ps p^ , which was not taken into account 
in the calculations. 
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FIG. 2. Osmotic pressure AP versus volume fraction (j> 
for latex A (a = 525 A). The symbols indicate experimental 
measurements from Benzing and Russel , for several initial 
concentrations n" of added electrolyte (NaCl) in the suspen- 
sion (mole/dm^): •, 0; O, 10"^ +, 10"^ A, 10""; □, 10"^ 
The arrow indicates the experimental point used to find the 
value of Za = 923. The solid curves represent our corre- 
sponding theoretical predictions for different values of initial 
concentration n° of added salt in the suspension (mole/dm'^): 
(a) 10"^; (b) 10"^ (c) 10""; (d) 10"^ We have not drawn 
the theoretical curve for n° = lO"'^ mole/dm^ since at this 
scale it is indistinguishable from the curve for n° = lO"'^ 
mole/dm"^. 
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FIG. 3. Osmotic pressure AP versus volume fraction <j) for 
latex Ji {a — 1830 A). The symbols indicate experimental 
measurements from Benzing and Russel for several ini- 
tial concentrations n^ of added electrolyte in the suspension, 
as identified in Fig. 2. The bar indicates the uncertainty in 
the experimental measurements. The arrow indicates the ex- 
perimental point used to find the value of Zb = 2212. The 
solid curves represent our corresponding theoretical predic- 
tions for different values of initial concentration of added 
salt in the suspension (mole/dm'^): (a) 10"^; (b) 10"^; (c) 
10"^ (d) 10"*. 

After the completion of this work, we became aware of 
the experimental measurements, performed by Reus et al. 



1^ , for a monodisperse charge-stabilized aqueous suspen- 
sion of bromopolystyrene particles. They also compared 
their measurements with the theoretical predictions using 
the Poisson-Boltzmann equation in a spherical WS cell 
1^, which is a calculation similar to ours of Section II, re- 
quiring the knowledge of the equilibrium concentration, 
PS, in the chamber containing pure electrolyte. Although 
this was explicitly provided by Reus et at, our formula- 
tion of Section III allows us to predict the osmotic pres- 
sure based on the a priori knowledge of the amount of 
salt that was being added to the system before the equi- 
librium is established. This is by far more practical than 
to wait for the establishment of the equilibrium and only 
then try to measure the salt concentrations in the two 
compartments and use this as an input into the Poisson- 
Boltzmann equation. Furthermore, Reus et al. provided 
the explicit value of the bare charge Z obtained using 
conductimetric titration and pH measurements. In Fig. 4 
we compare our theoretical predictions with the exper- 
imental data of Ref. A perfect fit is found without 
any need for adjustable parameter. 
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FIG. 4. Osmotic pressure AP versus volume fraction (jj 
for a monodisperse charge-stabilized aqueous suspension of 
bromopolystyrene particles of radius a — 510 A. The circles 
represent experimental measurements from Reus et al. 
obtained for an equilibrium concentration of salt in the pure 
electrolyte chamber ps = 10~^ mole/dm''. The bar indicates 
the uncertainty in the experimental measurements. The solid 
curve is our theoretical result obtained using the experimen- 
tally determined valence of the polyion, Z = 7500. Thus, in 
this case, there are no fit parameters. 

Finally, we should add that Benzing and Russel Q 
have also tried to estimate the value of the bare charge 
Z using the electrophoretic mobility at infinite dilution. 
This can be directly related to the surface potential of a 
colloidal particle by the Henry's law Given the sur- 
face potential, the bare charge of the particle at infinite 
dilution can be obtained through the direct use of the 
Debye-Hiickel theory. In Table I the bare charges found 
through this procedure, Z(est), are compared with the 
ones we found by fitting the data, Z(fit). We see that 



the values obtained through the use of Henry's law fall 
way below the values needed to fit the data. In view of 
our perfect agreement with the data of Reus et al. ^ , we 
are lead to the conclusion that there must be something 
wrong in trying to deduce the bare charge on the basis 
of electrophoretic mobility at infinite dilution. 



TABLE I. Comparison between the estimates for Z given 
by Benzing and Russel Z(est), and our values obtained 
through fit, Z(fit). 
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a (A) 


Z{est) 


Z(fit) 


A 


525 


174 


923 


B 


1830 


510 


2212 
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